Search for a 33.9 MeV/c 2 Neutral Particle in Pion Decay 
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The E815 (NuTeV) neutrino experiment has performed a search for a 33.9 MeV/c 2 weakly- 
interacting neutral particle produced in pion decay. Such a particle may be responsible for an 
anomaly in the timing distribution of neutrino interactions in the KARMEN experiment. E815 has 
searched for this particle's decays in an instrumented decay region; no evidence for this particle was 
found. The search is sensitive to pion branching ratios as low as 10 -13 . 

PACS numbers:14.80.-j, 12.60.-i, 13.20.Cz, 13.35.Hb 



The KARMEN collaboration at the ISIS spallation 
neutron facility at the Rutherford Appleton Laboratory 
uses a pulsed neutrino beam resulting from stopped pion 
and muon decays to study neutrino-nucleon interactions. 
Their experiment has reported an anomaly in the timing 
distribution of neutrino interactions from stopped muon 
decays jl). One possible explanation for the anomaly is 
an exotic pion decay, where a neutral weakly-interacting 
or sterile particle is produced and travels 17.7 m to the 
KARMEN detector with a velocity of 4.9 m//is. Upon 
reaching the KARMEN detector, the exotic particle de- 
cays to a partially electromagnetic state, such as e + e~v 
or 7za The e + e~v decay is strongly favored by recent 
KARMEN data ||. This slow moving exotic particle 
(hereafter denoted as Q°) would have a mass of 33.9 
MeV/c 2 , which is near the kinematic threshold for pion 
decay. 

The KARMEN experiment reports a signal curve for 
pion branching ratio B(ir — > /z + Q°) • B(Q° — > visible) 
versus lifetime. Their signal region extends as low as 
10~ 16 for a lifetime of 3.6 fxs. For branching ratios above 
this minimum, there exist two solutions to the KARMEN 
anomaly (at small and large lifetimes). Certain portions 
of the KARMEN signal have already been excluded. Ex- 
periments at PSI ||,|| have performed searches for this 
exotic particle by studying the momentum spectrum of 
muons and electrons produced by tt + decays in flight. 
PSI has excluded any exotic pion decays to muons with 
branching ratios above 2.1 x 1CP 8 at 90% C.L., and to 
electrons with branching ratios above 0.9 x 10~ 6 at 90% 
C.L. In addition, there exist astrophysical constraints on 



certain decay modes of the Q° which exclude lifetimes 
above 10 3 s ||. Despite the above limits, portions of the 
KARMEN allowed signal region remain to be addressed. 

The E815 (NuTeV) neutrino experiment at Fermilab 
has performed a direct search for the Q° decay by com- 
bining the capabilities of a high intensity neutrino beam 
with an instrumented decay region (the "decay chan- 
nel"). During the 1996-1997 fixed target run at Fermilab, 
NuTeV received 2.54 xlO 18 800 GeV protons striking a 
BeO target with the detector configured for this search. 
The secondary pions and kaons produced from the inter- 
action were subsequently sign-selected using a series of 
magnets and focused down a beamline at a 7.8 mrad an- 
gle from the primary proton beam direction. The pions 
and kaons could then decay in a 440 m pipe before hitting 
a beam dump. A total of (1.4 ± 0.1) x 10 15 pion decays 
and (3.6 ± 0.4) x 10 14 kaon decays occurred in the pipe. 
The neutral weakly-interacting decay products (neutri- 
nos and possibly Q s) traveled through approximately 
900 meters of earth berm shielding before arriving at the 
decay channel. 

The instrumented decay channel consisted of a series of 
helium bags, extending a total of 34 meters in length, in- 
terspersed with 3 m x 3 m multi-wire argon-ethane drift 
chambers. The drift chambers were designed to track 
charged particles from decays occurring within the he- 
lium. Upstream of the decay channel stood a 4.6 m x 4.6 
m array of scintillation plates, known as the veto wall, 
used to detect any charged particles entering from up- 
stream of the detector. Downstream of the decay chan- 
nel was the Lab E neutrino detector, which consisted 
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FIG. 1. Schematic of the NuTeV decay channel, including the veto wall, helium bags, drift chambers, and calorimeter. 



of a 690-ton iron-scintillator sampling calorimeter inter- 
spersed with drift chambers. The Lab E detector pro- 
vided triggering, energy measurement, and final particle 
identification for tracks entering from the decay chan- 
nel. Particles were identified by their penetration into 
the calorimeter: a muon produced a long track, a pion 
produced an elongated cluster of hits, and an electron or 
photon produced a compact cluster. More details of the 
decay channel and the calorimeter can be found elsewhere 

10- 

The decay channel was employed to search for the de- 
cay of the Q° . The experimental signature of this decay is 
a low-mass, low-transverse momentum electron-positron 
pair having a vertex within the fiducial volume of the 
decay channel. 

A series of analysis cuts isolated the Q° from back- 
ground. The cuts were divided into two categories: re- 
construction and kinematic. Reconstruction cuts isolated 
two track events occurring within the fiducial volume of 
the decay channel. We required that two charged tracks 
originated from a common vertex, with no additional 
tracks associated with the vertex. By removing events 
with activity in the veto wall, we ensured that no charged 
tracks entered from upstream of the decay channel. We 
also required that each track had a small slope (< 10 
mr) relative to the beam axis; and that, when projected 
upstream to the veto wall, it fell within 50" of the beam 
center. These cuts removed both cosmic rays and pho- 
tons from neutrino interactions in the upstream berm. 
Finally, we required that each track be identified as an 
electron based on the shower shape in the calorimeter. 
Because of the small opening angle of the two tracks, the 
two electron showers manifested themselves as a single 
merged electron-like cluster within the calorimeter. The 
efficiency for identifying such an event as an ee pair was 
estimated from Monte Carlo studies to be 90.1 %. The 
cluster energy was divided between the tracks based on 
a fit to the amount of multiple scattering each track un- 
derwent in the decay channel. 



Because of its low mass, the Q° possesses unique kine- 
matic features which can be used to distinguish it from 
potential background sources, such as photons and deep- 
inelastic neutrino interactions. We have used effective 
scaling variables to represent the kinematics of the re- 
constructed events. The effective scaling variables x c g 
and Wcff were calculated for each event using the fol- 
lowing assumptions: I) the event was a charged cur- 
rent neutrino interaction (ugN — > IN'X) and 2) the 
missing transverse momentum in the event was carried 
by an undetected final state nucleon. We have defined 

icff = 2m% is and w °fi = \j m l + 2rn p v v i s /c 2 - Ql is /c 2 , 
where Q v is is the visible reconstructed 4-momentum 
transfer, u V [ S is the reconstructed hadron (or electron) 
energy, and m p is the proton mass. Using the above defi- 
nitions, we required that all reconstructed events have 
x c ff < 0.001 and W eS > 2.5 GeV/c 2 . In addition, 
we required that the reconstructed transverse mass mx 
(rriT = \pr\ + yPrf m y > where ray is the invariant 
mass of the two charged tracks and pt is the momentum 
transverse to the beam axis) be less than 250 MeV/c 2 . 
Finally, we required that the total energy deposited by 
the e + e~ pair be greater than 15 GeV. The effect of these 
cuts when applied to signal and typical background kine- 
matic distributions can be seen in Figure | The total 
acceptance for Q° events in the fiducial region was 15.6%. 

The principal backgrounds originated from three main 
sources: neutrino interactions in the helium, neutrino in- 
teractions in the drift chambers, and neutral particles 
(mainly photons and kaons) from neutrino interactions 
in the berm and veto wall. Note that because the two 
electron tracks had a very small opening angle, the longi- 
tudinal vertex position resolution was quite poor (er w 7 
m). Thus, interactions in the chambers could not be 
removed using a vertex position cut. We used the Lund 
Monte Carlo program || to simulate the up deep-inelastic 
interactions in both the berm and the decay channel. 
Separate Monte Carlo programs were used to simulate 
hadronic resonance and diffractive single pion produc- 
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FIG. 2. Kinematic Monte Carlo distributions of x e s, visible 
energy, and transverse mass for neutrino interactions (solid) 
and Q° (dashed). All events shown have passed reconstruc- 
tion cuts only and are relatively normalized. 

tion MM. The GEANT Monte Carlo program (version 
3.21) pll| was used to simulate the decay channel and 
Lab E calorimeter. Based on our Monte Carlo study, we 
expected 0.06 ± 0.05 total background events within our 
signal region (See Table [j]). The uncertainty on this back- 
ground estimate is dominated by Monte Carlo statistics. 

A blind analysis was performed. The signal region was 
hidden while cuts were developed based on Monte Carlo 
studies of signal efficiencies and background rejection. 
Before examining the data near or in the signal region, 
we performed a series of studies to verify our background 
estimates. Using the Monte Carlo, we made predictions 
for the following three quantities: 1) the number of low 
energy (below 15 GeV) and high transverse mass (above 

TABLE I. Total Expected Background Events. Errors 
reflect Monte Carlo statistics. 



Source 


Rate 


Photons 


04+ 002 


Kaons 


< 


0.001 


Deep Inelastic Charged Current 


0.00 


± 0.04 


Deep Inelastic Neutral Current 


0.02 


± 0.02 


Cosmic rays 


< 


0.001 


Quasi-Elastic Charged Current 


0.000 


± 0.008 


Resonance Neutral Current 


0.000 


± 0.003 


Diffractive Pions 


± 


0.01 


Total 


0.06 ± 0.05 



TABLE II. Background Study Results. Uncertainties are 
systematic. 



Type of Event 


Events Predicted 


Events Seen 


High Transverse Mass 


2.0 ± 0.3 


1 


fin Events 


4.1 ± 0.6 


3 


Multiple Track Events 


13.7 ± 1.8 
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FIG. 3. Kinematic distributions of multiplicity, energy, and 
transverse mass for data (crosses) and background Monte 
Carlo (dashed) multi-track events. Monte Carlo is absolutely 
normalized. 

500 MeV/c 2 ) background events; 2) the number of fnr 
events; and 3) the number of multi-track events occur- 
ring within our decay channel. The results of these stud- 
ies (shown in Table |l|) demonstrate good agreement be- 
tween data and the Monte Carlo predictions. In the case 
of multi-track events, where a larger sample of events 
was available, there was also good agreement for various 
kinematic distributions (see Figure y) . 

The above acceptance estimate was based on a heavy 
neutrino model where the Q° decays to eev |Q . We also 
considered additional decay models where the Q° decays 
to and 771/. The acceptances under these decay sce- 
narios were (0.5 ± 0.1)% and (1.1 ± 0.1)%, respectively. 
The low efficiencies are due to the requirement that two 
tracks be reconstructed; this was only possible if a photon 
converted in the low-mass decay region before entering 
the calorimeter. 

Systematic errors for this result were dominated by 
uncertainties on the number of pion decays in the pipe 
(6.8%) and the overall normalization. The sensitivity 
normalization was taken from a measurement of the num- 
ber of neutrino interactions in the decay channel using 
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FIG. 4. Branching ratio versus lifetime plot for the Karmen 
signal and the exclusion regions ( 90% C.L.) from the NuTeV 
and PSI experiments. Systematic errors (except for decay 
model) have been included. 



very loose cuts. This number was (10.0 ± 4.3)% below a 
prediction normalized to the number of neutrino interac- 
tions in the calorimeter. Considering this disagreement 
as a systematic error on the normalization, we have calcu- 
lated a total systematic error of 12.1% on the sensitivity. 

Upon analyzing the signal region, we found no events 
which passed the selection criteria. The probability of 
seeing zero events from an expected background of 0.06 ± 
0.05 is 94%. We thus present an upper limit, shown in 
Fig. |, on B(ir -> n + Q°) ■ B(Q° -> visible) || NuTeV is 
also sensitive to ir — > e + Q°; the limit is shown in Fig. ||. 

This result excludes a region of parameter space which 
extends as low as four orders of magnitude below current 
limits on the short lifetime solution to the KARMEN 
anomaly. An experiment with significantly more pion 
decays will be necessary to confirm or rule out the longer 
lifetime and lower branching ratio regions. 

This research was supported by the U.S. Department of 



* A simple interpretation of the Q° is that the particle is a 
sterile heavy neutrino which mixes with the muon neutrino. 
This requires a branching ratio B(n — > fj, + Q°) k6x 1(T 8 , 
which has already been ruled out by the PSI limit Still 
viable is a model where the heavy neutrino is produced by 
a smaller mixing with z/ M and decays primarily via a larger 
mixing with v T , such that KQVm)! • l(Q°K)l ~2x 1CT 6 . 
Within this model, we set a limit on this product " 
|(Q°|^ T )| < 1.4 x 10~ 3 at 90% C.L. 
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FIG. 5. Branching ratio versus lifetime plot for the NuTeV 
limit for ty — > e + Q°. Systematic errors (except for decay 
model) have been included. 
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